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We consider the forward-backward pion charge asymmetry for the e"'"e~ 7r"'"7r~7 process. In 
addition to the Bremsstrahlung and the double resonance contributions, we use four different models 
to describe the final state radiation (Kaon Loop Model, Resonance Chiral Perturbation Theory, 
Unitarized Chiral Perturbation Theory and Linear Sigma Model), perform a Monte Carlo calculation 
»' ' and compare our results with experimental data. Our results indicate that none of these models 

, yield an appropriate description of the asymmetry in the whole energy range reported by the KLOE 

collaboration and, contrary to expectations, the invariant mass region between 400 and 700 MeV 
{^JQ' may be the region of major interest since the models we consider are unable to reproduce the data 

'~i , in this region while the same models (except for Resonance Chiral Perturbation Theory) reproduce 

. the overall characteristics of the data in the /o(980) region. 
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I. INTRODUCTION 



The nature of low mass scalar mesons nonet is a long-standing puzzle. The 4> radiative decays are expected to provide 
(~| ' information about the /in_(980) and an (980) scalar mesons. Unfortunately data reported by the KLOE collaboration 
'— ' on the (f> decays to /o7 [ill and ao7 H - with tt^tt^j and 7r°7?7 final states respectively- together with results on the 
^ ' process 4> 7r+7r"7 including the /07 as intermediate state [3| are not conclusive. In the latter work, results on the 
^ , forward-backward asymmetry as a function of the 7r"'"7r~ invariant mass are presented. The asymmetry is sensitive to 
' the mechanisms involved in the final state radiation [3| and has been proposed as an appropriate observable to test 
the importance of the different mechanisms at work in this reaction and the different proposals for their description. 
T-H It also provides information on the pion form factor [5,]. 

(N 

QQ The asymmetry requires a non vanishing interference between initial (ISR) and final (FSR) state radiation, the latter 
being strongly model dependent |6|]. The invariant amplitude for the e^e~ —^ 7r"'"7r~7 process can be parameterized in 
0^ . terms of three independent Lorentz structures and thus the model dependence in FSR can be included in three scalar 
■ functions The final state radiation has been calculated in different models. The simplest approximation has 

J> , been named scalar QED (sQED) [6] and it actually includes the p contributions to the pion form factor. In that work 
the contributions of intermediate scalars (/o(980) and cr) are also considered using a point-like 0/o7 interaction, in the 
so called "no-structure" model. Later on, final state radiation was calculated within Resonance Chiral Perturbation 
Theory (RxPT) [3 at tree level [i,0,@- In particular, in d] sub-leading intermediate vector mesons contributions 
like e~^e~ ^ — » p~tt~^ —t tt'^tt~j, named double resonance contributions, were incorporated. 



The aim of this paper is to work out the predictions for e+e^ —>■ tt^tt^j at the (j) resonance at one loop level for 
four alternative models, namely RxPT, Unitarized Chiral Perturbation Theory (UxPT) [l^ (containing actually a 
resumation of loops) , Linear Sigma Model (LSM) [il|,d2i] and the so-called "kaon-loop" model [3]. In each case we 
add the sQED * VMD model and the double resonance contribution, both proposed in Q. We report the forward- 
backward pion charge asymmetry which is calculated using Monte Carlo methods and compare our results with KLOE 
data. 



The paper is organized as follows: Section II includes the general formalism to describe the e^e^ —^ tt^tt^j process. 
In section III we derive the scalar functions fi that characterize the RxPT, LSM, UxPT and the kaon loop (KL) 
models, including Bremsstrahlung process and the double resonance contribution. In section IV we present the 
numerical results after the Monte Carlo calculation and compare them with experimental results. Finally, conclusions 
are given in section V. 



2 



II. GENERAL FORMALISM 



We are interested in the process 

e" (pi) e+ (P2) ^ 7r+ tt" (p-) 7 (fc, e) • (1) 

For completeness, in order to introduce our notation and conventions, in this section we include the basic equations 
used to describe the process. To this end, we follow the formalism developed in Ref. J|. The invariant amplitude AI 
includes the initial state radiation Misr, and final state radiation AIpsR, M = Misr + Mfsr, with 



MisR = -^L^^^ell^F, , (2) 



, ^r<- (3) 

where F.^ (g^) denotes the pion electromagnetic form factor, t^, is the photon polarization vector and the tensor Mp'^ 
describes the photon radiation from the final state. The lepton currents are given by 



^2 tl 



xwsi(pi), (4) 



= eus^ {-P2) lt,Us, {pi) . (5) 

The electron and positron spinors are u^^ (pi) and (— P2) respectively. In terms of the external particles' four- 
momenta, the following variables are introduced Q = pi + p2,q = p+ + P--,l = p+ — P- and five independent Lorentz 
scalars are defined 

s = Q2 = 2pi ■p2, 

ti = {Pi - kf -ml = -2pi ■ k, 

t2 = {P2 - kf = -2p2 ■ k, (6) 

Ui=l-pi,U2 = l-p2, 

where the electron mass, me, has been neglected. The differential cross section is 



with p+ = (i?-|_,p+) = (i?_,p_), k ^ [uj = \k\ , k) and is the invariant amplitude squared, averaged over 

initial lepton polarizations [24] . The most general form of the FSR tensor Mp^ is 7j 

Mr - h^" + f2rr + /3<^ («) 

where the r/"^ are three independent gauge invariant tensors which are dictated by parity, charge conjugation, crossing 
symmetry and gauge invariance 

= k^Q-' - g^-'k ■ Q, 

= A; • ? (Z^Q'' - .g^'^fc • + (fc^fc -l-Fk-Q), (9) 
Tg^" = g2 (g^^-'fe . / - fcA'/'') + (/^fc ■Q-Q'^k-l). 

The scalar functions /,; = /i (Q^, k ■ Q,k ■ are either even (/i,2) or odd (/a) under the change of sign of the argument 
k-l. Our first task will be to determine these scalar functions fi for RxPT, U xPT, LSM and the KL model including 
Bremsstrahlung process and the double resonance contribution f§| . Before we consider the models, we define the charge 
asymmetry. 

The pair of pions produced in ([1]) differ in charge conjugation, depending if the photon is emitted from the initial or 
from the final state, while the former is odd under charge conjugation the latter is even. So, any interference between 
the two amplitudes is odd under charge conjugation and gives rise to a charge asymmetry. The forward-backward 
charge asymmetry is defined as 

_ 7V(g^+ > 90°) - jV(g^+ < 90°) 

N{e^+ > 90°) + N{e^+ < 90° ^ ' 

where 0^+ is the 7r+ polar angle, which is measured with respect to the incident electron momentum. It should be 
clear that the asymmetry depends strongly on the experimental conditions, in particular on the cutoff polar angle 
and the minimal photon energy that can be measured. 
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FIG. 1: Feynman diagrams for the 'ycpmr vertex function in RxPT 

III. FSR MODELS 

A. Bremsstrahlung process and double resonance contribution 

Before discussing the (j) decay models, we first consider the Bremsstrahlung and the double resonance contribution. 
The Bremsstrahlung process was calculated in Q usingthe sQED * VMD model and RxPT. The functions fi for 
this contribution are given by equations (II) to_(20) in [9|. The double resonance contribution e+e ~> ch p tt^ — > 
7r+7r~7, is formulated in the context of VMD [ij] and the fi terms can be found in equations (26) to (28) in We 
used these expressions in our numerical calculation. We now consider specific models. 

B. The RxPT and UxPT models 

The amplitude for the e^e^ — > 7r+7r^7 process at the (f> peak involves the ^ (p-K'^ it" vertex function with all particles 
on-shell. On the other hand, within the analysis of e+e^ — > ^tttt reaction, this vertex function was calculated in the 
context of RxPT and U xPT ^15] for an off-shell photon. The starting point was the one loop calculation of the j(j>Tnr 
vertex function [l^l which turns out to be finite. The diagrams relevant to our analysis are shown in Fig. ([p. These 
diagrams include kaons in the loops, thus they involve the off-shell KK — tttt amplitude. It was shown in that, 
to leading order in the chiral expansion, the contribution of diagrams d,e,f,g cancels the off-shell contributions of 
the KK ~ TTTT amplitude, entirely contained in diagrams a, b, c, h, so that the calculation reduces to evaluate diagrams 
a,b,c,h with the KK — mr amplitude on- shell. This procedure yields the RxPT result for the 707r"'"7r~ vertex 
function and we would expect it to reproduce experimental results at low dipion invariant mass. However, due to the 
appearance of the widely discussed light scalar resonance (the a meson), this expansion brakes down in the scalar 
channel even at the dipion threshold. 

The scalar poles can be generated unitarizing the leading order meson-meson scattering amplitudes for definite isospin. 
Following 16], the unitarized K~^K~ — tt+tt" scattering amplitude is calculated projecting onto the zero spin and 
isospin channel the leading order on-shell meson-meson amplitude and performing a coupled channel analysis involving 
iterations of all intermediate states in the s channel. 

As far as the calculation of the e+e~ — > 7r"''7r~7 amplitude is concerned, the scalar poles are incorporated replacing 
the leading order on-shell K^K^ — tt+tt^ amplitude by the K^K^ — tt+tt" unitarized amplitude. For details of 
the calculation we refer the interested reader to [IS], here we just quote the result for 7r+7r~7. The resulting 
amplitude for 0(Q, r]""") tt+ {p+) tt" (p_) 7 (k, e^) in UxPT is 
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-iM = 



V3 2V2TT^mlP V3 



Gv [lp{Q -kg^,,- Qf,K)^ Qa 



(11) 



where = {p+ and t\^^ denotes the unitarized isoscalar scalar KK — rnr amplitude. The factor 2/\/3 is 

required to single out the 7r~'~7r~ contribution in the isoscalar tttt channel. The function gKio^) is given by 



9K{qn = -1 + log^ 



o{q^ 



(12) 



1 



3^. Note that the particular form of this function involves a regularization scheme as well as 



with a[q^) - 

a subtraction point. We use dimensional regularization and the value /i = l.2GeV which reproduces the /o peak at 
9S0MeV in the squared meson- meson amplitudes in the scalar channel 16]. The loop integral is given by 



rab 



/(^) 




(13) 



z > 



z < 



Results for RxPT are obtained from Eq. pT|) by replacing t^j^^ by the leading order on-shell interaction Vktt = 
~\f%q^ I^Lp . Using the propagator for a vector meson in the tensor formalism we obtain the amplitude for 7*((3, [>) 
7r7r7(fc, v) via the exchange of the vector meson as 



iM 



le 



■FvV2 



Q2 - M' + iV^M^ 



where 



/ 
J 



Gv 



''Ktt Tab 



Gv ~ 



Fv\ mj. 



9K{q^)- 



In terms of this vertex function we can identify the final state radiation invariant tensor Mp'^ (see Eq. [8]) 

M^"- = -*e2(ATr + /2rr + = -^e'M^- 



with 



1 Fv 1 (Q^ 



, ah 
yip 



■Gvl 



G\ 



2 



9K{q^) , 



/l 

/3 = 0. 

Notice that Eq. fTS)) contains a term with the combination Gv — ■ This combination is small and it vanishes in 
the context of Vector Meson Dominance [l^l ■ For the sake of completeness we will keep this term and study its effect 
on the charge asymmetry. 



(14) 

(15) 
(16) 

(17) 

(18) 

(19) 
(20) 
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C. The phenomenological KL model 

In this model the process under consideration proceeds through the chain 

e- (pi) e+ (p2) ^ ^ /o (g) 7 (fc, e) - 7r+ {p+) (j>-)l (k, e) • (21) 
The corresponding amplitude is 

■ 2 
— ?p 

= Aevip2) i^u (pi) (Q^'fc'' - Q • kg^'") e,, (22) 



where we have defined 



with 



A = '-f^TT^^lfF, {s) (q^) e»^-, (23) 



and f^s, (70, and stand for the /ott+tt^, (pK^K^ , /qK^K^ and — 7 couphngs respectively (for details concerning 
the precise definition of these quantities we refer the reader to the appendix of Ref. flS||). The angle is the phase 
of elastic background and it must be included with the kaon loop factor This phase becomes relevant for 

the interference between ISR and FSR amplitudes and thus can be measured for charged pions in this case 
5b = hyjq^ - 4m2 with h = 757GeV [13, The kaon loop function if is given in Finally, the scalar 

functions for this model are obtained by comparing ((22|) with ([3]) and ([8]), in this way we get 

/i = A 

h = 0, (25) 
/3 = 0. 



D. The Linear Sigma Model 



The calculation in this approach is similar to the KL model, the difference arising from the treatment of the scalars. 
For the neutral pion case, the amplitude has already been derived in [l^ using the improved chiral loop approach. 
Thus we can obtain the amplitude we are interested in just by making the following replacement in the KL amplitude 



9s9fPf {q') - A {K+K- ^ 
where the amplitude for the meson scattering is given by 



= V2A(K+K- 



^0 0^ 



LaM 



(26) 



AlK+K- 



9 9 



LaM 



K 



K 



Da (g2) 



-C0S ccps 



V2t 



S(ps 



'if 



m 



S(^s ( s^s + V2c^s 



(27) 



with Ds (q^) = q^ — m'g + imsTs, Jk = 1-22/^, (c0s, s^s) = (cos sin^s) and 05 is the scalar mixing angle in the 
strange- non-strange basis [l^- Therefore, the scalar functions for this model can be obtained from (I^S)) by making 
the replacement indicated in the previous paragraph. 



IV. NUMERICAL RESULTS 



In order to obtain the numerical results we implemented a Monte Carlo code where the experimental conditions of the 
KLOE collaboration are included. Thus, for the 7r+ polar angle - defined respect to the electron beam - we considered 
the range 45° < 9t^+ < 135°. As far as the photon is concerned we took 45° < 0^ < 135° and used a cut at 10 
MeV, i.e. we assumed > 10 MeV [3|. In order to compare with the experimental data for the forward-backward 
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500 550 600 

Invariant mass (MeV) 



FIG. 2: Comparison of models predictions with KLOE data in region I (400-700 MeV). A scalar mixing angle (j>s ~ —5° is 
used in the LSM calculations. The predictions of the four models are in strong disagreement with data. 



asymmetry it is convenient to consider three complementary regions in the dipion invariant mass. We observe that in 
the 400- 700 MeV the data decreases slowly while in the 700 - 900 MeV range the asymmetry grows up. Above 900 
MeV, at 930 MeV, the asymmetry decreases until it reaches 950 MeV where it starts growing again and reaches the 
maximum in the /o(980) region. Hereafter we will refer to these as regions I, II and III respectively. In the following 
paragraphs we compare our numerical calculation with the experimental data. 

Calculations in UxPT and RxPT involve parameters that have already been fixed from meson phenomenology. 
We use the following values in the numerics: Gy = 53 MeV, Fy = 154 MeV, U = 93 MeV and ^ = 1.2 GeV [T&|. As 
to the KL model, a summary of the values we use for the parameters of this model is given in Table [H including the 
Qs and gf values reported in and the Particle Data Group [2l|. The complete definitions of these constants are 
described in p^ . 



Parameter 


Value 


Reference 


TUf (MeV) 


980 


[21] 


Ff (MeV) 


70 


[21] (Averaged) 


g's (GeV^) 


3.61 


[20] 


2 
94. 


19.56 


[13, 21] 


gj (GeV^) 


7.78 


[20] 




179.14 


[13, 21] 


9l 


35.95 


[13, 21] 


fl 


24.66 


[13, 21] 



TABLE L Constants appearing in the KL model and the numerical values used in this work. 

Besides the intrinsic parameters of the scalar mesons, the LSM involves the scalar mixing angle. We use the following 
values: fUf = 980 MeV, T/ = 70 MeV, while for the the sigma meson we took the values reported in [l^l rria = 528 
MeV and To- = 414 MeV and for the scalar mixing angle we consider three values: 4>s = —3° (LSM3), —5° (LSM5) 
and -7° (LSM7). 

In the following lines we briefly describe comparison of KLOE data with the model predictions for the forward 
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FIG. 3: Comparison of the forward backward asymmetry as predicted by RxPT and UxPT with KLOE data. The points 
labeled UxPTw exclude contributions from the Gv — ^ term in Eq. (|18p. Region II (left) and region III (right) 



backward asymmetry (Regions: I 400 - 700 MeV, II 700 - 900 MeV and III above 900 MeV): 

• Fig. ([2]) contains our results for the four models in region I. We observe that in this region the results of all four 
schemes are in strong disagreement with the experimental data. 



• Figure ([3]) presents the UxPT and RxPT results for regions II and III. The difference between UxPT and 
UxPTw results is that the latter exclude contributions from the Gy — ^ term in Eq. (|18p (VMD limit |23]). 
We observe that in region II both models provide a satisfactory description of the behavior of the asymmetry. 
In region III U xPT predictions are in very good agreement with data. The effect of the Gy — ^ term in Eq. 
(fT8)) is small but necessary to achieve a good description of data. In this region RxPT does not reproduce the 
data. 

• Kaon loop model predictions for regions II and III are shown in figure (j4|). In region II the results are system- 
atically above the data, while in region III the predictions of this model follow the trend of the data. 

• LSM predictions for regions II and III are shown in figure In region II good agreement with data is obtained 
for the three values of the mixing angle. The asymmetry in region III is highly sensitive to the value of the 
scalar mixing angle and a very good description of data is obtained for 05 = —5°. This value is close to the one 
reported in |19j where the /o propagator is approximated by a Breit-Wigner. 

• Figure ^ contains a summary of the best predictions of the four models we consider in this work. 



A. Phase dependence of the asymmetry. 

The strong disagreement between the models and the KLOE data for low dipion invariant mass motivated us to 
consider the effect of a relative phase. We explore the effect of a constant relative phase between sQED * VMD 
plus double resonance model contributions and the models considered here. We included a constant phase e*" in the 
KL, LSM and UxPT amplitudes. The calculation was performed for different values of a and the results are shown 
in figures ([7]) and Results are labeled with the corresponding phase, e.g. LSMA5 corresponds to the Linear 
Sigma Model with a phase a = 45°. An strong dependence on the relative phase is observed. In the case of UxPT a 
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FIG. 4: Comparison of KLOE data with the results for the forward-backward asymmetry in the kaon loop model . Region 11 
(left) and region III (right). 




700 720 740 760 780 800 820 840 860 880 900 

Invariant mass (MeV) 



< 

CQ 0,24 



—\ ' r- 



■ 


KLOE 


• 


LSM3 


▲ 


LSM5 


T 


LSM7 



920 940 960 980 

Invariant mass (MeV) 



1000 1020 



FIG. 5: Comparison of KLOE data with results for the forward-backward asymmetry as predicted by the LSM using the scalar 
mixing angles (f)s = -3° (LSM3), -5° (LSM5) and -7° (LSM7). Region II (left) and region III (right). 



relative phase close to a = 135° yields a reasonable description of region I. However, such a phase destroys the good 
description of the data in Region III. The inclusion of a phase in the KL and LSM does not improve the description 
of data. We conclude that a constant relative phase does not provide a good description of data in the whole invariant 
mass range reported by the KLOE collaboration for the models considered here. 
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FIG. 6; Comparison of KLOE data with predictions of the four models. Region II (left) and region III (right). 
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FIG. 7: Phase effect on the forward-backward asymmetry within the KL (left) and LSM with <f>s ~ —5° (right). 



V. SUMMARY AND CONCLUSIONS 



We used Unitary Chiral Perturbation Theory, Resonance Chiral Perturbation Theory, Linear Sigma Model and 
the Kaon Loop model to describe the final state radiation of the e+e" 7r+7r~7 process at the loop level. We 
implemented the models in a Monte Carlo code in order to obtain the forward-backward asymmetry from these models 
and the conventional tree level mechanisms {sQED * VMD plus double resonance model). Detailed comparison with 
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FIG. 8: Phase effect on the forward-backward asymmetry for UxPT. 



experimental data is reported. The following conclusions can be drawn form the comparison of the models with the 
KLOE data: 

• Predictions of the four models for the asymmetry in the low invariant mass region (400 - 700 MeV) are in strong 
disagreement with the KLOE data. 

• In region II (between 700 and 900 MeV) the four models yield an appropriate description of KLOE results. 

• The trend of the data in region III (above 900 MeV) is well reproduced by all models except for RxPT. 

• The asymmetry in region III is highly sensitive to the scalar mixing angle in the Linear Sigma Model and a 
value 0s = —5° is favored by data. 

• None of the models studied here yields an appropriate description of the asymmetry in the whole range of energy 
reported by the KLOE collaboration. 

• The inclusion of a relative phase between tree level contributions and the loop level contributions studied here 
is not helpful. In the UxPT case a phase close to a = 135° improves the prediction in the low invariant mass 
but this phase destroys the good description of data in the /o(980) region. 

• Contrary to expectations, the invariant mass region between 400 and 700 MeV may be the region of major 
interest. In fact, models are unable to reproduce the data in this region while the same models correctly 
reproduce the overall characteristics of the data in the /o(980) region. A possible interpretation is that the 
/o(980) pole is correctly described in these models while the joint description of poles in different regions 
requires a more careful analysis. 
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